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WHAT IS CLAIMED: 

1 . A method for predicting pressure independent permeation 
flux and target molecule 3deld in a permeate resulting from crossflow membrane 
filtration of particles in a poly-disperse suspension, said method comprising: 

determining particle size distribution of the poly-disperse 

suspension; 

determining equivalent spherical radii of the particles; 
determining viscosity of the suspension; 

detenniningiiiaximum back-transport velocity (Mi) for all particles; 

estimating maximum aggregate packing volume fraction for 
all particles at a wall of the filtration membrane from geometric considerations; 

selecting the particle that gives a minimnm permeation flux at a 
given filtration membrane shear rate, wherein the selected particle has a radius 

determining a predicted permeation flux (J)\ 

determining packing density at a membrane wall for each 
particle size (aj for / ^ i) at the predicted permeation flux; 

determining interstitial packing density i^untersdce) of particles in 
the suspension which are the smallest; 

determining minimum pore diameter (2rmmimuiii) based on the 
packing density of each particle; and 

estimating yield of a target species in the filtration permeate by 
calculating observed sieving coefficient (So) for the target species, thereby 
predicting permeation flux and target molecule yield of the poly-disperse 
suspension during crossflow filtration. 

2. The method according to claim 1, wherein said determining 
viscosity of the suspension is carried out by using a modified Einstein- 
Smoluchowski equation: 7/70 = 1+ 2,5(i%,+ ki^fe^ where rj is bulk fluid viscosity 
(kg/m.s) of the suspension, 70 is bulk fluid viscosity of the suspension without 
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solute (kg/m-s), h is particle shape factor (-), and <fe is particle volxime fraction in 
the bulk suspension (-). 

3- The method according to claim 1, wherein said determining 
viscosity of the suspension is carried out by experiment 

4. The method according to claim 1, wherein said determining 
maximum back-transport velocity (wi ) comprises: 

calculating Brownian difiEiision (Jb) for all particles, where 

calculating inertial lift (Ji) for all particles, where Ji— 

0.036pa^//Tj; 

calculating shear induced diffusion (Js) for all particles where Js = 
0.07S(a^/Ly^^jiix{^^i,% wherein y is wall shear rate (s"^), kis Boltzmaim constant 
(J/mol K), ris temperature (K), tj is bulk fluid viscosity Otg/m.s), at is radius of 
species i(m), L is tube length (m), ^ is particle volume fraction at the filtration 
membrane (-), ^ is the particle volume fraction in the bulk suspension (-), and p 
is particle density (kg/m^); and 

selecting /max for each particle, wherein Jhiax = ui whereby 
maximum back-transport for each particle is determined. 

5. The method according to claim 1, wherein said estimating 
maximum aggregate packing voliraie fraction(<Jji^ at the membrane wall for a 
suspension comprises: 

determining particle size (ai) of species (i) in the suspension; 

detemiining if the size ratio of the particles is > 10, such that a^-i > 
lOajforall^z;; and 

calculating the maximum aggregate packing volume fraction( jJm),) 
by i^Afii = + (1- ?W«-7), where ^M=^is set to 0.64 when the size ratio of the 
particles is > 10, such that a^+i > 10a, forall a,-. 
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6. The method according to claim 5, wherein the suspension 
comprises 3 particle sizes and wherein ai >10a2 > lOOaa, said method further 
comprising: 

calculating ^a/ = ^ + (1- <4„) + 0.74[1 - + (1- ^)}], 
wherein ^ is set to 0.64. 

7. The method according to claim 1, wherein said estimating 
maximum aggregate packing volume fraction at the membrane wall for a 
suspension comprising two particles, such that ai >10 ai. is carried out by 
calculating ^m = ^ + 0.74 (1- ^n^), where is set to 0.64. 

8. The method according to claim 1, wherein said estimating 
TnaTTiTnuTTi aggregate packing volume fraction at the membrane wall 
comprises: 

calculating a maximum radius ratio of all particles; 
determining if said maximum radius ratio is < 10; and 
setting (^as 0.68, where said maximum radius ratio is < 10. 

9. The method according to claim 1, wherein said selecting the 
particle that gives the minirmim permeation flux (J) comprises: 

calculating Brownian diffusion (Jb) for all particles, where 
JB=0Al4ip^I^/7fa^Lf^]ni^J(pi;); 

calculating inertial Ufl (Ji) for all particles, where Ji= 

calculating shear induced difiusion (/s) for all particles, where Js = 
0.078iaVLy^]in{^^b\ wherein y is wall shear rate (s"^), Kis Boltzmann constant 
(J/mol K), ris temperature (K), tj is bulk fluid viscosity (kg/m.s), a,- is radius of 
species z (m), L is tube length (m), ^ is particle volume fraction at the membrane 
wall (-), is the particle volume fraction in the bulk suspension (-), and p is 
particle density (kg/m^); 

determining a Jmax value for each particle; and 
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selecting a Jmax value from among all Jmox values that is tiie 
lowest, th^eby selecting the minimum pemieation flux (J). 

10. The method according to claim 1 , wherein said determining 
packing density at the membrane wall {jpy^j) for all particles at the predicted 
permeation flux (ay for j ^ i) comprises: 

back-calculating the value of ^ such that gives the predicted 
permeation flux (J) of selected particle (aO using the equation for back-transport 
that establishes maximum back transport for each particle (ay fory=0, wherein the 
equation is either Jb = OA\Aiyf^T^lTj^(^Lf^\n{(lKj^b) or 

Js = Q.01Ka'^ILf^'Mi<f^lh\ or Ji= ^mepc^-^lri, where y is wall shear rate (s"^), 
K\s Boltzmaim constant (J/mol K), Tis temperature (K), rj is bulk fluid viscosity 
(kg/m.s), a, is radius of species z(m), L is tube length (m), ^ is particle volume 
fraction at the membrane wall (-), <fe is the particle volume fraction in the bulk 
suspension (-), and pis particle density (kg/m^). 

1 1 . The method according to claim 10, wherein said 
determining packing density ftirfher comprises: 

determining if the predicted penneation flux is established by 
inertial lifl: for one particle type; 

determining if uji >10J; and 

setting =0^ when one particle type is established by inertial lift 
(7i)andwyi>10J: 

12. The method according to claim 10, wherein said 
determining packing density fiirther comprises: 

determining if the predicted permeation flux is established by 
inertial lift (J{) for one particle type; 

determining if uji <10J; and 

determining packing density (^J) by ^yi = ^w^^j when uji 
<10Jand one particle type is established by inertial lift. 
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1 3 . The method according to claim 1 0, whttein said 
determining packing density further comprises: 

determining if pameation flux is established by inertial lift {J\ ) for 

more than one particle type; 

determining if wji < 10/ for the particles; and 

determining packing density by = $W - S when uji < lOJ 

and peraieation flux is established by inertial lift for more than one particle type, 

14. The method according to claim 1 0, wherein said 
detennining packing density ftirther comprises: 

detCTmining if permeation flux is established by Ji for more than 
one particle type (jll, jI2,. . jin); and 

determining packing density at the membrane wall by 
<*Mdii = ^Af- 2 ^y, wherein : = Uya : ^ji2 

where j ^gll or jI2 and mjh, Uji2< lOJi when permeation flux is established by /i 
for more than one particle type. 

1 5. The method according to claim 1 , wherein said determining 
interstitial packing density i^unterstice) of the smallest particle is carried out by 

^Unterstice = ^mcorrected/i\ " ^ ^corrected^^ whcrcin 4K^icorrected = $W [ 

where ^ is flie particle volume fraction at the membrane wall (-) for particle t 

16. The method according to claim 1, wherein said determining 
minimum pore diameter (2rmiiiimuni) is carried out using 

2rnuniiiium= «i{V2[4(4/3)7i/^//„^^^^ -2}, where a is radius of 
species i (m) and rnrimmumis a minimum equivalent cake void radius for all cake 
types (m). 

1 7. The method according to claim 1 , wherein said estimating 
yield of a target species comprises: 

calculating observed sieving coeflScient (iS^^, where *So = SJUl- 
Sa)exp('J/k) + Sa), wherein actual sieving coefficient Sa is obtained from 



wo 2004/016334 



PCTAJS2003/025230 



-61 - 

Sa = (iS'ocexp(jPem )y(Sa: + exp(Pem ) -1), woU Peclet ntimber, Pem is obtained from 
Peja = (JSJD)(Sa: I s<f> whcTe Jis permeation flux (m/s), 5n is taken as the 
side of the &ce centered cube of the particles of radius ax that forms the 
controlling cake for transmission^ ^ = a = ai [(4(4/3)7i;)/ifert/errtfce]^'^9 D is 
molecular difiEiision coefficient (m^/s), intrinsic sieving coefficient is obtained 
from &c = (1-A,)^[2 - (1- X)^] exp(-0,7146A?), A, = rjr^ where rs is solute radius 
(m) and is a minimTm equivalait cake void radius for all cake types (m), ^ is 
equilibrium partition coefficient between membrane pore and suspension (-), s is 
cake/membrane porosity (-), is hindrance factor for diffiisive transport (-), and 
k is mass transfer coefficient (m/s). 

18. The method according to claim 1, wherein crossflow- 
filtration is carried out in a diafiltration mode, and the yield of the target species is 
estimated after iVd diavolumes as Yield =1- exp(- JVdiSoaverage) where iSoaveragc is 
average observed sieving coefficient during diafiltration (-), where ~ SJUl' 
Sa)o^p(rJ/k) + Sa ), where actual sieving coefficient Sa is obtained from Sa = 
(iScexpCPem )y(Saz + expCP^m ) —1), where J is permeation flux (m/s), wall Peclet 
number, Pe^ay is obtained from Pem, = (J^D)(Soz ls(j>Ki)^ where ^ is taken as tiie 
side of the face centered cube of the particles of radius ai that forms the 
controlling cake for transmission, where <Sn = ^ = «i M^M^I'^y^l ^imterstux^^^ ^ D is 
molecular diffusion coefficient (m^/s), intrinsic sieving coefficient Ac is obtained 
from iSoc = (l-Xf\2 - (1- Xf^ exp(-0.7146A,^), X = rjrram, where rs is solute radius 
(m) and rmin is a minimum equivalent cake void radius for all cake types (m), ^ is 
equilibrium partition coefficient between membrane pore and suspension (-), s is 
cake/membrane porosity (-), is hindrance factor for diffixsive transport (-), and 
k is mass transfer coefficient (m/s). 

19. The method according to claim 1 further comprising: 
re-calculating packing density for all particle sizes if packing 

constraiQts are not satisfied based on initial determination of packing densities of 
the particles at the wall. 
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20. The method according to claim 19 further comprising: 
correcting packing density using (^icorrected = <kf iiMI^ M\ 
reevaluating J for the particle selected as having the minimum 

permeation flux based on t^orrected- [ $^/]; and 

reevaluating maximum back-transport velocity (wi)* 

2 1 . The method according to claim 20 further comprising: 
repeating the st^s of claim 17 until a desired packing constraint is 

met 

22. The method according to claim 1 further comprising: 
refining the yield of the target species. 

23. The method according to claim 22, wherein said refining 
the yield comprises: 

determining whether the suspension has a low, intermediate, or 
high operating shear rate leading to different yield regimes, wherein a suspension 
at a low op^ting shear rate leads to an 5^ ^ 0.75 corresponding to a yield > 0.95, 
an intermediate operating shear rate leads to 0 < iSo < 0.75 corresponding to yield 
fi-om 0 to 95%, or a higji operating shear rate leads to an iSi, s 0, wherein 
Sa/iil- iSa)exp(-J/t) + Sa)^ wherciu actual sieving coefficient Sa is obtained firom Sa 
= (Soc&xp(Pem ))/(5x + expCP^m ) -1), wall Peclet number, Pejn is obtained fi-om 
Pcm = (JSm/D)(Soc /s^Kd), where J is permeation flux (m/s), ^ is taken as the side 
of the face centered cube of the particles of radius ai that forms the controlling 
cake for transmission, 4i = ^ = «i [i'^(^/^)^y^inie,^tice\^^ D is molecular diffusion 
coefficient (rr^/s), intrinsic sieving coefficient Saz is obtained firom Sec = (1-A,)^[2 - 
(1- Xf] exp(-0.7146;i^), X rs/rnun, where rg is solute radius (m) and r„^ is a 
tniiiimum equivalent cake void radius for all cake types (m), ^ is equilibrium 
partition coefficient between membrane pore and suspension (-), s is 
cake/membrane porosity (-), Kd is hindrance factor for diffusive transport (-), and 
k is mass transfer coefficient (m/s). 
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24. The method according to claim 23, wherein an intermediate 
operating shear rate is determined as leading to 0 < iSq < 0.75, said method further 
comprising: 

calculating stagnant film flux (J) equation for non-retentive 
membranes wherein 7= k In [(<|)wi - <|>penneatei) - <l>penneatd)] = In - 
5o)], wherein (<|)wi » <|>t)enneatdi); and 

correcting So by replacing J- solvent permeation flux ( m/s) with 
the stagnant film flux {J) equation for non-retentive membranes in the equation for 
observing sieving coefficient, iSi,, where So = SJ((1' Sa)exp(rJ/lc) + Sa ). 

25. The method according to claim 1 further comprising: 
constructing a plot of the predicted permeation flux and yield 

versus wall shear rate, thereby predicting permeation flux and target molecule 
yield of the poly-disperse suspension during microfiltration. 

26. The method according to claim 1, wherein filtration is 
selected firom the group consisting of microfiltration and ultrafiltration. 

27. The method according to claim 1 , wherein filtration is 
carried out with a SLtec selected firom the group consisting of a flat sheet filter, 
hollow-fiber filter, and a helical filter. 

28. The method according to claim 1, wherein the suspension is 
selected from the group consisting of streams firom biomedical and bio-processing 
industries, waste water, surface water, environmental pollutants, industrial waste 
streams, and industrial feed streams. 

29. The method according to claim 28, wherein the suspension 
is a stream firom biomedical and bio-processing industries selected firom the group 
consisting of proteins, cells, nucleic acids, colloids, milk, and suspended particles. 
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30. A method for determining packing density of particles of a 
poly-disperse suspension at a membrane wall, said method comprising: 

providing a predicted permeation flux (J); 

determining packing density for all particle sizes at the predicted 
permeation flux; and 

detemuning interstitial packing density {^untenticd of particles in 
the suspension which are smallest, thereby detemiining packing density at tiie 
membrane wall of particles of the poly-disperse suspension. 

3 1 . The method according to claim 30, wherein said 
determining packing density at the membrane wall for all other particles at 
the predicted permeation flux (ay for j ^ 0 comprises: 

back-calculating the value of (j>wj such that <pwj gives the predicted 
permeation flux (J) of selected particle (aQ, using the equation for back-transport 
that establishes maximum back transport for each particle {aj for y=0> wherein the 
equation is either Jb = 0.1 14(;7S£^I^/77^a^X)^^ln(^/^) or 

Js = 0.078(a'*/i)^^;in((^/<^b), or /i= OSBepc^^lrj, where y is wall shear rate (s"^), 
Kis Boltzmann constant (//mol K)^ Tis temperature (K), tj is bulk fluid viscosity 
(kg/m,s), ads radius of species z(m), L is tube length (m), ^ is particle volume 
fraction at the membrane wall (-), ^ is the particle volume fraction in the bulk 
suspension (-), and pis particle density (kg/m^). 

32. The method according to claim 3 1 , wherein said 
determining packing density further comprises: 

detemiining if tihie predicted permeation flux is estabUshed by 
inertial lift (J{) for one particle type; 

determining if uji >10J; and 

setting j^wy = 0, when one particle type is established by inertial lift 

33 . The method according to claim 3 1 , wherein said 
determining packing density further comprises: 
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detennimng if the predicted permeation flux is established by 
inertial lift (J\) for one particle type; 

determining if uji <\QJ; and 

determining packing density (^) by ^ji = $W - 2 when uji 
<\QJ and one particle type is estabhshed by inertial lift, 

34. The metibiod according to claim 3 1 , wherdn said 
determining packing density fiirther comprises: 

determining if pemieation flux is established by inertial lift {Ji) for 

more than one particle type; 

determining if mji < lOJfor the particles; and 

deteraiining packing density by =^ (pM-^ <fKvj when wji <\0J 

and permeation flux is established by inertial lift for more than one particle type. 

35. The method according to claim 3 1 , wherein said 
determining packing density ftirther comprises: 

detennining if permeation flux is established by Ji for more than 
one particle type (jll » jI2,. . .jlw); and 

determining packing density at the membrane wall by 

+ = - 2 wherem : = <feju Wji2 : jfejE «ju, 
where j 9fcjll or jI2 and wjn, Wji2< 1(W when permeation flux is estabhshed by Ji 
for more than one particle type. 

36. The method according to claim 30, wherein said 
detamioing interstitial packing density i^onterstice) of the smallest particle is 
carried out by 

(1 - E (fK^jcorrected^^ Whcrcin Corrected - ^ [ (M^^ 

where ^ is tihie particle volume fraction at the membrane wall (-) for particle /. 

37. The method according to claim 3 1 fiirther comprising: 
re-calculating packing density for all particle sizes and 
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deterraining if packing constraints are not satisfied based on initial 
determination of packing densities of the particles at the wall. 

38. The method according to claim 37 further comprising: 
correcting packing density by nsmg ^iconected ^ ^m\. {^rL $^/]; 
reevaluating Jfor the particle selected as having the minimum 

permeation flux based on ^comcted^ ^ [ (^f)/2 j^f]; and 

reevaluating maximum back-transport velocity (Mi). 

39. The method according to claim 30, wherem filtration is 
selected firom the group consisting of microfiltration and ultrafiltration. 

40. A method for predicting pressure independent permeation 
flux for crossflow membrane filtration of a poly-disperse suspension, said method 
comprising: 

determining viscosity of the suspension; 

determining maximum back-transport velocity (wi ) for aU particles; 

estimating maximum aggregate packing volume firaction {^m) for 
all particles at a wall of the filtration membrane from geometric considerations; 

selecting the particle that gives a minimum permeation flux at a 
given filtration membrane shear rate, wherein the selected particle has a radius 

determining a predicted permeation flux (J); and 

determining packing density (^y) at the membrane wall for each 

particle size (a, for j ^ i) at the predicted permeation flux, thereby predicting 

pressure independent permeation flux for the suspension. 

41 . The method according to claim 40 further comprising: 
re-calculating packing density for all particle sizes if packing 

constraints are not satisfied based on initial determination of packing densities at 
tiie wall. 
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42. The method according to claim 41 further comprising: 
correcting packing density using (jhmcorrected = [ 

reevaluating / for the particle selected as having the minimum 
permeation flux based on ^wicomcted^ $W [ and 

reevaluating maximum back-transport velocity (ui). 

43 • The method according to claim 40, whorein said 
determining viscosity of the suspension is carried out by using a modified 
Einstein-Smoluchowski equation: 7/770 = H- 2.5 + ki^^, where tj is bulk fluid 
viscosity (kg/m.s) of the suspension, 70 is bulk fluid viscosity of the suspension 
without solute (kg/m.s), ki is particle shape factor (-), and ^ is particle volume 
fraction in the bulk suspension. 

44. The method according to claim 40, wherein said 
deterniining viscosity of the suspension is carried out by experiment 

45. The method according to claim 40, wherein said 
determining maximum back-transport velocity (wi) comprises: 

calculating Brownian diffusion (/b) for all particles, where 

Jb =0.114(y3^7^/7^a^Ly^^hir<(>H/j*6); 

calculating inertial lift (/i) for all particles, where Ji= 

0.036pa^//7; 

calcxdating shear induced diffusion ( Js) for all particles, where Js = 
0.078(a'*/L)^^;in(^Jw/$^b), and wherein y is wall shear rate (s**^), isris Boltzmann 
constant (J7mol K), Tis temperature (K), rj is bulk fluid viscosity (kg/m.s), a,- is 
radius of species /(m), L is tube length (m), ^ is particle volume fraction at the 
membrane wall (-), is the particle volume fraction in the bulk suspension (-), 
and pis particle density (kg/m^); and 

selecting ./max for each particle, wherein /max = Ui ,thereby 
determining maximum back-transport for each particle. 
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46. The method according to claim 40, wherein said estimating 
inf^yiTTniTri aggregate packing volume fi:action(^Ai) 3t the membrane wall for a 
suspension comprises: 

determining particle size (ai) of species (0 in the suspension; 
determining if the size ratio of the particles is > 10, such that an-i > 
10a; for all a,-; and 

calculating the maximum aggregate packing volume fi:action((i^ 

by 

= + (1- where set to 0.64, when the size 

ratio the particles is >10, such that an-i > lOa,- for all a,-. 

47. The method according to claim 40, wherein the suspension 
comprises 3 particle sizes and wherein ai >10a2 > lOOaa, said method further 
comprising: 

calculating ^ = ^, + <2>m (1- M + 0-74[l - {^^ + (1- ^i)}], 
wherein is the yns^Tcmnm packing volume fraction for monodisperse spheres set 
to 0.64. . 

48. The method according to claim 40, wherein said estimating 
maximum aggregate packing volume fraction {^j^ at the membrane wall 
comprises: 

calculating a maximimi radius ratio of all particles; 
determining if said maximum radius ratio is < 10; and 
setting 0.68, where said maximimi radius ratio is < 10. 

49. The method according to claim 40, wherein said estimating 
maximum aggregate packing volume fraction ( j*Ai) at the membrane wall for a 
suspension comprising two particles, such that ai >10 ai, is carried out by 
calculating = + 0.74 (1- where ^ is set to 0.64. 

50. The method according to claim 40, wherein said selecting 
the particle that gives a rniiiimum permeation flux (J) comprises: 
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calculating Brownian diffusion (Js) for all particles, where 
calculating inertial lift (/{) for all particles, where Ji= 

0.036/wV/77; 

calculating shear induced difiusion (Js) for all particles, where Js = 
0,078(a^/Ly^jia(^^, wherein y is wall shear rate (s"^), isris Boltzmann constant 
(J7mol K), ris temperature (K), tj is bulk fluid viscosity (kg/te.s), ai is radius of 
species i(m), L is tube length (m), is particle volume fraction at the membrane 
wall (-), is the particle volume fraction in the bulk suspension (-), and p is 
particle density (kg/m^); 

determining a Jmax value for each particle; and 
selecting a Jmax value from among all Jmax values that is the 
lowest thereby selecting the minimum permeation flux (J). 

5 1 . The method according to claim 40, wherein said 
determining packing density at the membrane wall {^J) for all particles at the 
predicted permeation flux (ay for j ^ i) comprises: 

back-calculating the value of such that gives the predicted 
permeation flux (J) of selected particle (aj), using the equation for back-transport 
that establishes maximvim back transport for each particle (aj for y =0? wherein the 
equation is either Jb = 0.1 lA^j^I^/Tfa^Ly^hii^/^b) or 

Js = 0.07S(aVLy^jiai^/M^ or Ji= 0.036pa^//7j, where y is wall shear rate (s"^), 
^ris Boltzmann constant (J/mol K), Tis temperature (EC), 7j is bulk fluid viscosity 
(kg/m.s), a,- is radius of species /(m), L is tube lengfii (m), ^ is particle voliime 
fraction at the membrane wall (-), is the particle volume fraction in the bulk 
suspension (-), and p is particle density (kg/m^). 

52. The method according to claim 5 1 , wherein said 
determining packing density further comprises: 

determining if the predicted permeation flux is established by 
inertial Uft (70 for one particle type; 

determining if uji ^1 OJ; and 
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settmg = 0, when one particle type is established by inertial lift 
and uji >\QJi. 

53 . The method according to claim 5 1 , wherein said 
determining paddng density fiarther comprises: 

detenniiiing if the predicted peroieation flux is established by 
inertial lift (Ji) for one particle type; 

detemiining if uji <\0J; and 

determining packing density (^) by = h^-^ ^ when uji 
<1(X/ and one particle type is established by inertial lift. 

54. The method according to claim 51, wherein said 
determining packing density further comprises: 

determining if permeation flux is established by inertial Uft (Ji) for 

more than one particle type; 

determining if Wji < lOJ^for the particles; and 

detemiining packing density by $jwji = - E when t/ji < 10/ 

and permeation flux is established by inertial lift for more than one particle type. 

55. The method according to claim 51 , wherein said 
determining packing density ftirther comprises: 

determining if permeation flux is estabUshed by Ji for more than 

one particle type (jll, jI2,. . jin); and 

determining packing density at the membrane wall by 

<*Mdn + ^vidi2 = $W - 2 wherein <l>y^vi : = ^Ajii Wji2 : <^&ji2 Wjii, 

where j ^gll or jI2 and Wju, mje < 10^ when permeation flux is estabUshed by Ji 

for more than one particle type. 

56. The method according to claim 40, wherein filtration is 
selected ft^om the group consisting of microfiltration and vdtrafiltration. 
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57. A method for calculating yield of a target molecule in a 
permeate for a poly-disperse suspension during crossflow mCTobrane filtratioxi, 
said method comprising: 

determining minimum pore diameter (2rmiiiiimiin) based on the 
packing density of each particle and 

estimating yield of a target species in the filtration permeate by 
calculating observed sieving coefficient (So) for the target species. 

58. The method according to claim 57, further comprising: 
refining the yield and pressure independent permeation flux. 

59. The method according to claim 57, wherein said refining 
the yield comprises: 

determining whether tbe suspension has a low, intermediate, or 
high operating shear rate leading to different yield regimes, wherein a suspension 
at a low operating shear rate leads to an iS^ > 0.75 corresponding to a yield > 0.95, 
an intermediate operating shear rate leads to 0 < < 0.75 corresponding to yield 
from 0 to 95%, or a high operating shear rate leads to an So= 0, wherein 5^,= 
SJiiU Sa)exp(rJ^) + 5^), wherein actual sieving coefficient Sa is obtained firom Sa 
= (iSQcexp(Peni)y(«S'oc+ exp(Pem) -1), wall Peclet number, Pem is obtained firom 
Peja = (J^D)(Soc /s^Kd)^ where J is permeation flux (m/s), <Sn is taken as the side 
of the face centered cube of the particles of radius ai that forms the controlling 
cake for transmission, wh^e ^ = a = ai [(4(4/3)n)/^iintersticef^, D is molecular 
diffusion coefficient (m^/s), intrinsic sieving coefficient 5^ is obtained from iSoc = 
(1-A,)^[2 - (1- Xf] exp(-0.7146X.^), X = rjrrr^m. where rs is solute radius (m) and rmin 
is a minimum equivalent cake void radius for all cake types (m), ^ is equiUbrium 
partition coefficient between membrane pore and suspension (-), s is 
cake/membrane porosity (-), Kd is hindrance factor for diffusive transport (-), and 
it is mass transfer coefficient (m/s). 
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60. The method according to claim 59, wherein an intemiediate 
operating shear rate is determined as leading to 0 < iSq < 0.75, said method further 
comprising: 

calculating stagnant film flux (J) equation for non-retentive 
membranes wherein J = k In [(<|)wi - <|j^cnneatei) K^i - <l)penneatci)] = ^ In l^b^^ - 
iSo)], wherein » <|>^enneatei); and 

correcting Sq by replacing J— solvent permeation flux ( m/s) with 
the stagnant film flux (7) equation for non-retentive membranes in the equation for 
observing sieving coefficient, iSo, where So = SJifl- Sa)oxp('J/k) + Sa ). 

6 1 . The method according to claim 57, wherein determining 
minimum pore diameter (2rminimuin) is carried out using 

2rniimmum = «i{^2[4(4/3)7i/^vviinteretice]^^ -2}, wbere a is the radius of 
species i (m) and rmimmumis a minimum equivalent cake void radius for all cake 
types (m). 

62. The method according to claim 57, wherein said estimating 
yield of a target species comprises: 

calculating observed sieving coefScient (SoX where 
So = SJai' Sa)expi-J^) + Sa ), wherein actual sieving coefficient Sa is obtained 
from Sa — (SccGxpiJPeta ))/('S'oc + exp(Pem ) —1), wall Peclet number, Pem is obtained 
from Pem = (JStiJD)(Soc /s^Ka)y where J" is permeation flux (m/s), is taken as 
the side of the face centered cube of the particles of radius ai that forms the 
controUiag cake for transmission, where <5bi = a = [(4(4/3)K)/^iintentice]^^^^ D is 
molecular diffusion coefficient (w^/s% intrinsic sieving coefficient Soc is obtained 
from Soc = (l'^f[2 - (1- Xf] exp(-0.7146A,^), X = rs/rmm, where rs is solute radius 
(m) and is a minimum equivalent cake void radius for all cake types (m), ^ is 
equilibrium partition coefficient between membrane pore and suspension (-), s is 
cake/membrane porosity (-), is hindrance factor for diffusive transport (-), and 
^is mass transfer coeffident (m/s). 
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63. The method according to claim 57, wherein crossflow 
filtration is carried out in a diafiltration mode, said the yield of the target species 
after Nd diavolumes is estimated by Yield =1- exp(- iVd5i)average)j where iSoavaage is 
average observed sieving coefficient during diafiltration (-), where So = SJUl' 
Sa)G}qp{-J/k) + Sa ), where actual sieving coefficient Sa is obtained firom Sa = 
(iSocexpCPem ))/(5'oc + expiPcm ) -1), wall Peclet number, Pcja, is obtained firom Peja 
= (J^D)(Soz where Jis pemieation flux (m/s), <Sii is taken as the side of 
the face centered cube of the particles of radius oi that forms the controlling cake 
for transmission, where 4i = « = «i [(4(4/3)%)/ ^untersticef^, D is molecular 
difiusion coefficient (po^/s), intrinsic sieving coefficient Soc is obtained firom Sec = 
il'Xfll - (1- Xf] Qxp(-0J146X\ X = rs/rmin, where rs is solute radius (m) and rmin 
is a minimum equivalent cake void radius for all cake types (m), and ^ is 
equilibrium partition coefficient between membrane pore and suspension (-), € is 
cake/membrane porosity (-), is hindrance factor for diflEusive transport (-), and 
^is mass transfer coefficient (m/s), 

64. The method according to claim 57, wherein filtration is 
selected firom tiie group consisting of microfiltration and ultrafUtration. 

65- The method according to claim 57, wherein filtration is 
carried out with a filter selected from the group consisting of a flat sheet filter, 
hollow-fiber filter, and a helical filter. 

66. The method according to claim 57, wherein the suspension 
is selected from the group consisting of streams from biomedical and bio- 
processing industries, waste water, surface water, environmental pollutants, 
industrial waste streams, and industrial feed streams. 

67. The method according to claim 66, wherein the suspension 
is a stream from biomedical and bio-processing industries selected from the group 
consisting of proteins, cells, nucleic acids, colloids, milk, and suspended particles. 
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68. A method for designing a orossflow membrane fQtration 
system for a poly-disperse suspension, said method comprising: 

selecting a poly-disperse suspension; 

applying the method according to claim 1 to predict pressure 
independent permeation flux and target molecule yield in a permeate for the 
selected poly-disperse suspension; and 

optimizing conditions for filtration based on the prediction of 
pemieation flux and target molecule yield to design a filtration system for the 
selected poly-disperse suspension. 

69. The method according to claim 68, wherein filtration is 
selected firom the groiip consisting of microfiltration and iiltrafiltration. 

70. The method according to claim 68, wherein filtration is 
carried out with a filter selected firom the group consisting of a flat sheet filter, 
hoUow-fiber filter, and a helical filter. 

71 . A method of selecting operating conditions of a crossflow 
filtration system for poly-disperse suspensions, said method comprising: 

applying the method of claim 1 to a crossflow filtration system for 
a selected poly-disperse suspension to determine a limiting pressure independent 
permeation flux for a given shear rate and expected yield of a target species for 
the selected system conditions and 

selecting the operating conditions of the system using the 
determined limiting pressure independent permeation flux for a given shear rate to 
obtain an optimal balance between permeation flux and yield of a target species. 

72. The method according to claim 7 1 , wherein filtration is 
selected firom the group consisting of microfiltration and ultrafiltration. 



73. The method according to claim 71 , wherein filtration is 
carried out with a filter selected firom the group consisting of a flat sheet filter, 
hollow-fiber filter, and a helical filter. 
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74. The method according to claim 71, wherein the suspension 
is selected from the group consisting of waste water, surface water, environmental 
pollutants, industrial waste streams, and industrial feed streams. 

75. A method of modeling a process for filtration of a poly- 
disperse suspension comprising: 

applying the method according to claim 1 for a poly-disperse 
suspension using a computer-generated program to model a process for filtration 
of the poly-disperse suspension. 

76. The method according to claim 75, wherein filtration is 
selected from the group consisting of microfiltration and ultrafiltration. 

77. The method according to claim 75, wherein filtration is 
carried out with a filter selected from the group consisting of a flat sheet filter, 
hollow-fiber filter, and a helical filter. 

78. The method according to claim 75, wherein the suspension 
is selected from the group consisting of streams from biomedical and bio- 
processing industries, waste water, surface water, enviroimiental pollutants, 
industrial waste streams, and industrial feed streams. 

79. The method according to claim 78, wherein the suspension 
is a stream from biomedical and bio-processing industries selected from the group 
consisting of proteins, cells, nucleic acids, colloids, milk, and suspended particles. 



